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Abstract 

The properties of warm symmetric and asymmetric nuclear matter 
are investigated in the frame of the Thomas- Fermi approximation us- 
ing a recent modern parametrization of the effective nucleon-nucleon 
interaction of Myers and Swiatecki. Special attention is paid to the 
liquid-gas phase transition, which is of special interest in modern nu- 
clear physics. We have determined the critical temperature, critical 
density and the so-called flash temperature. Furthermore the equation 
of state for cold neutron star matter was calculated. 

PACS: numbers: 21.65.+1, 21.10.Dr, 21.60.-n, 25.75-q 



*This investigation is dedicated to Georg Siifimann at the occasion of his seventieth 
birthday. 
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I Introduction 



An essential but also a very complex and complicated problem of modern 
physics is the structure of matter under extreme conditions of temperature 
and/or density. The equation of state of hot dense matter is a key ingredient 
for many branches of modern physics. To mention in this respect are, for 
instance, the understanding of the physics of the early universe, supernova 
explosions and protoneutron stars , respectively. A further important exam- 
ple are high-energy heavy-ion collisions, which provide, at present, the only 
laboratory means to study the behaviour of hot and dense matter. In the 
literature one can roughly distinguish between two theoretical attempts to 
investigate the equation of state, namly the non-relativistic and relativistic 
approaches. The majority of the relativistic treatments are performed in the 
framework of the relativistic Hartree approximation and only a few investiga- 
tions use the relativistic Hartree-Fock or relativistic Briickner-Hartree-Fock 
approximation (see, for instance, Refs. 0-0). The majority of the theo- 
retical treatments utilize the non-relativistic scheme, using either effective 
density dependent interactions or the Bruckner approach 0-0. 

In this investigation we are going to extend a recent modern Thomas- 



Fermi approach of Myers and Swiatecki to the equation of state at 

finite temperatures (see also Ref. jn]]). Furthermore we will briefly discuss 
the question of neutron star matter and the equation of state for higher 
asymmetries, respectively. 

The paper is organized as follows: First, we recapitulate fore the sake of 
completeness in Section [Ij the basic structure of the Thomas-Fermi model 
|12| , |l3f for density and momentum dependent interactions. The extension 



of the model to finite temperatures is given in Section III. Section IV deals 



with the results of the model and the last Section is a summary of the results 
with conclusions. 



II Thomas- Fermi model at T = 

We shall not recapitulate the basic assumptions and limitations of the Thomas- 
Fermi model in nuclear physics, since they are described in greater detail in 
the pioneering work of Myers and Swiatecki [JT1| . In the modern version they 
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use an effective interaction of the following structure (r denotes the isospin): 
2T FjT fr 12 \ , , 

Vl2r = 9 (HI) 



Po V a 

i (lTf)a _i (lTC) (,M 2 




with the following definition: 

f$ = l(/i+f$), (n.2) 

here, pi, p 2 denote the (proton or neutron) density at the locations T\ and 
r 2 , respectively. The density at equilibrium is p = ( t^o) -1 ' an d the Fermi 
momentum and the Fermi kinetic energies are given by [m := |(m n + m p )\. 

Pf = » (|t 2 Po) * , (113) 

T F = |, (II.4) 

^ - ft < IL5 ' 
By means of the parameters £ and £ the interaction distinguishes the forces 
between nucleons of the same kind (upper sign) and of different kind (lower 
sign). With the parameter a one adjusts mainly the binding properties of 
nuclear matter. The repulsion is described by the momentum dependent 
term oc fip\ 2 - Since this repulsion turns out to be too strong for higher 
relative momenta, one corrects this deficit by the term oc TIP12 1 — 1 - 111 order 
to obtain better values for higher asymmetries, the authors have chosen £ 7^ C 
in contrast to former structures of Seyler and Blanchard fTSfl. The term 



proportional to <r(2p/ po) 2 ^ 3 takes care for a better agreement with the nuclear 
optical potential (for more details, see Refs. [|71-|]TT]). Fore our task the model 
contains six free parameters (in the expression for the energy the parameters 
(3 and a enter only in the combination B = (3 + |cr) and they are given in 
Table I. For the space dependence the standard Yukawa form is employed, 
i.e.: 

r 12 1 exp(-^f) 

9{ - ] = w^^' (IL6) 
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The energy per nucleoli in the Thomas-Fermi model is given by: 
E 1 2 



u 



N p (2tt^) 3 
with the single-particle potential: 



+ (II ' 7) 



Vr(pi) = -T^Ni ( / ^P2^12r + J d 3 p 2 V 12 ) , (II.8) 



{2nTi) 

or more explicitly: 

2 2 



(2tt^) 3 po 



T f ,t dp2\ai — Pi[ — + 7j | r - (Ti — 

^ \ VPf/ IP12I VPo 

+T F / d 3 p 2 f a u - A* ( — ] + 7«7^t - o"« ( — 



(II.9) 



Here, «i, f3\, 71, crj and a u , /3 U , j n , a u are defined as: ot\ = 0.5(1 — £)a; 
A = 0.5(1 - C)/3; 7i = 0.5(1 - 07; *i = 0.5(1 - CV; « u = 0.5(1 + £)a; 
Ai = 0.5(1 + C)/3; 7u = 0.5(1 + C)7 and a u = 0.5(1 + Qa. 

The results of this model for cold matter are described in more detail 
in Refs. [f7|-|]ll| and agree rather well with the values of the semiempirical 



droplet mass formula. Interesting is also the equation of state (EOS) for cold 
neutron star matter, which is characterized by the generalized /3-equilibrium 
between nucleons and leptons, i.e. electrons and muons (for more details, see 
Ref. fll6|). One has to obbey in this case the additional constraint for the 



chemical potentials, i.e.: 

£ip = fin - Pe] £ie = /V, (11.10) 

and charge neutrality between the protons and the leptons: 

P P = Pe-+P^-- (11.11) 
The EOS for this case is contained in Figs. 1 and 2. 
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Ill Thomas- Fermi model at T > 



The generalization of the Thomas-Fermi model for finite temperatures is 
formally more or less staightforward, since one has mainly to incorporate 
the Fermi distribution functions into the scheme for zero temperature. One 
obtains then for the energy per nucleon: 



u 



+00 , 

E / d 3 Pl (^ + -V T ( Pl ))f T (p,), (IIL12) 



3 

T 



with the temperature dependent single particle potential: 
2 2Xp T 



V T (pi) = 



(27Tft) 3 po 



*M*-a(S ,+ ^-«(S>w 



— oc 

2 2T F 



(27T^) 3 po 

+00 / / \ 2 



(111.13) 

and the density is now given by: 



2 



E / d 3 Pl f T (pi). (HI.14) 



(2tt^) ; 

The Fermi distribution function, / r , is defined as (k^ = 1): 

/r(Pi) = (l + exp (^(e T ( Pl ) - p' T )))) 1 , (111.15) 
where e T denotes the single-particle energy: 

2 

^(pi) = /^ + K(pi). (111.16) 
2m T 
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With the single particle energy one can rewrite the energy per nucleon as: 



u 



i + °° 1/2 \ 

^¥- P vl dp M^ + ^ M ) jM - (iili7) 



Due to the momentum dependence of the single particle potential one can 
also introduce an effective mass via: 

m . = fl^MV 1 . (IIL i8) 



p dp J 



p=0 



For finite temperatures one needs the free energy in order to obtain the 
pressure and the chemical potentials. The free energy per baryon is given 
by: 

f = u-Ts, (111.19) 

with the entropy per baryon, s: 



E / rf 3 p 1 (/.(Pi)ln/.(pi) + (l-/.(pi))ln(l-/ r (p 1 )) N 

(111.20) 

The expressions for the chemical potentials and the pressure are: 

l>r=\p(-S-) +11/ (III.21; 



P-r,T 



(fi' T in Eq. |III.15| is different to fi T , since the interaction is density dependent, 
see Appendix A of Ref. [0]) and 

p= 'S*(£L- (IIL22) 



More explicit expressions are given in Ref. JTT 



IV Results 

First we show in Figs. 1 and 2 for the sake of completeness the EOS for cold 
nuclear matter with different asymmetries [5 := (p n — p p )/p] and neutron 
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star matter. In Figs. 3 and 4 the EOS for symmetric nuclear matter (5 = 0) 
and pure neutron matter (5 = 1) is displayed for different temperatures. For 
small densities, below p = 0.0165 fm~ 3 , one obtains as expected, the be- 
haviour of a free Fermi gas with a linear temperature dependence, since the 
nucleon-nucleon force has a small range. For increasing density the EOS ex- 
hibits a quadratic temperature dependence, since in the temperature range 
smaller than the Fermi temperature a quadratic temperature dependence 
u(T) = u(T = 0) + a v T 2 holds (see Refs. 0, |i~7|). The chemical potential 
p, the pressure, the free energy per baryon and the entropy per baryon for 
symmetric matter are exhibited in Figs. 5, 6, 7 and 8 for different temper- 
atures. The values of the entropy per baryon are in good agreement with 



experimental results (see Ref. p0|), this is important, since the entropy per 



baryon is also a nessesary ingredient for the calculation of protoneutron stars 



(see e.g. Refs. p], 0). 

Due to the non-monotonic behaviour of the isotherms one can infer that 
these systems undergo a liquid-gas phase transition. Such a transition, which 
may occur in warm matter produced in heavy-ion collisions, is one of the most 
interesting problems in modern nuclear physics. Several groups are investi- 
gating this problem by studying multifragmentation, where the coexistence 
of nuclear fragments and free nucleons could be considered cLS cL SI ernal for 



the transition phase (see, for instance, Refs. [[H], p2|). The critical tem- 
perature can be obtained by a Maxwell construction. The situation can be 
either viewed in the P-p-diagram or in the p-P-diagram, respectively. Both 
views are given for symmetric nuclear matter in Figs. 9 and 10, respectively. 
Gibbs condition of phase equilibrium demands that pressure, temperature 
and chemical potential are equal in the two phases. 

As can be seen from the EOSs the nuclear matter incompressibility de- 
creases with increasing temperature. As a consequence one should exspect 
a limiting temperature, at wich a self-bound system can still exist in hy- 
drostatic equilibrium (P = 0). Beyond this temperature, the so-called flash 



temperature Tfl, the unbound system starts expanding p3| . The isothermal 
incompressibility is defined as: 



2«7 . „_8/ 



and is zero at the point where the minimum pressure reaches the zero pres- 
sure line (here at a density of p = 0.0973 fm -3 ). We obtain for the flash 
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temperature the value 15.61 MeV which is, as expected, smaller than the 
critical temperature T c , because at T c the pressure is already positive (see 
Fig. 6). If one increases the asymmetry of the system the critical temperature 
and the critical density decrease. These dependencies are depicted in Figs. 11 
and 12. For asymmetries larger than 0.853 the liquid-gas phase coexistence 
disappears in our model. 



V Conclusions 

We have calculated the EOS for cold and warm nuclear matter for different 
asymmetries using the approach of Myers and Swiatecki. Almost all non- 
relativistic calculations give critical temperatures in the range of 14 - 22 
MeV and critical densities of approximately 1/3 of the saturation density 0, 
which is in accordance with our results. Investigations for finite temperatures 
are relatively scare and the model seems to provide a simple theoretical tool 
for investigating the nuclear EOS in such ranges, which seems to be of greater 
interest, since it becomes now possible by means of the heavy-ion physics to 
study the behaviour of matter different from the saturation point, which is an 
essential ingredient in the theory of neutron stars and supernova explosions 
and also in the interpolation of heavy-ion collisions . 

In our opinion, it seems worthwhile, due to its simplizity and exellent 
agreement with nuclear data, to use this model with possible improvements 
in the further non-relativistic investigations of the EOS. Also the properties 
of neutron stars seem to be in accordance with this approach WE . 
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Table captions 

TABLE I: Parametrization of the interaction 
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TABLE I 



Parameter x Value Value x\ Value x. 



a 


0.59542 fm 


a 


3.60928 


P 


0.37597 


1 


0.21329 


a 


1.33677 


B 


1.48995 




0.44003 


c 


0.59778 



1.01054 2.59874 

0.07561 0.30036 

0.04289 0.17040 

0.26884 1.06793 

0.29964 1.19030 
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Figure captions 

Fig. 1. Energy per nucleoli versus density for nuclear matter for different 
asymmetries and neutron star matter in /3-equilibrium. The nuclear 
matter parameter are: po — 0.16545 fm~ 3 ; u = 16.527 MeV; incom- 
pressibility K = 301.27 MeV; asymmetry parameter J = 31.375 MeV. 
The maximum star mass for a non-rotating neutron star is approxi- 
mately 2.0 solar masses and the star radii are between 10 and 12 km. 

Fig. 2. Pressure as function of density for cold nuclear matter for different 
asymmetries and neutron star matter. 

Fig. 3. Energy per nucleon for symmetric nuclear matter (5 = 0) versus 
density at different temperatures. 

Fig. 4. Energy per nucleon neutron matter (5 = 1) versus density at differ- 
ent temperatures. 

Fig. 5. Chemical potential as function of density for symmetric nuclear mat- 
ter at different temperatures. 

Fig. 6. Pressure versus density for symmetric nuclear matter at different 
temperatures. 

Fig. 7. Free energy per baryon versus density for symmetric nuclear matter 
at different temperatures. 

Fig. 8. Entropy per baryon versus density for symmetric nuclear matter at 
different temperatures. The entropy behaviour agrees with the experi- 
mental situation (cf. with Fig. 6 of Ref. ||, see also Ref. p0|). 

Fig. 9. Liquid-gas phase transition in the P-p diagram for symmetric mat- 
ter: Shown are the isoterms for T = 0, 10, 15.61, 20.8 and 30 MeV. 
Above the dashed-dotted curve the matter exists only in one phase, 
either as gas (g) between the dashed-dotted line and the isotherm with 
the critical temperature T c = 20.8 MeV or as liquid (f) above T c = 20.8 
MeV. Between the dashed-dotted curve and the dashed curve one finds 
metastable phases, again either as a gas (g') or a liquid (P). Below the 
dashed curve the state is not stable. Shown is also the pressure at the 
flash temperature (15.61 MeV). 
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Fig. 10. Chemical potential as function of the pressure for symmetric nu- 
clear matter. The crossing point on each isotherm is the condition of 
phase equilibrium. 

Fig. 11. Asymmetry dependence of the critical temperature T c . 

Fig. 12. Asymmetry dependence of the critical density p c . 
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FIGURE 1 
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FIGURE 3 




-20.0 1 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 

p[fm 3 ] 



FIGURE 4 
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FIGURE 5 
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FIGURE 7 
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FIGURE 9 




FIGURE 10 
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